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Abstract 



Timing of Geminga gamma-ray pulsar is done using data of COS B and EGRET. It is shown, that 
errors in angular coordinates of sources similar to Geminga strongly influence a determination of i>, so 
that at angular precision less then 10~^ arc sec determination of the value of i> by means of criterion 
gives error more then 100%. Attempt have been done to improve coordinates of a gamma-ray pulsar 
using timing analysis. In addition to search of i/, v and D, a technique is first developed permitting 
search of two angular coordinates, absolute speed value and direction of a proper motion. In that way 
timing of gamma-ray pulsars gives amount of information compatible with radiopulsars, but using of 
' real data gives much poorer precision. In gamma-ray sources with rare pulses the periodicity criteria are 

I quite different from the ones in radio region. Data on the coordinates and proper motion of Geminga, 

obtained from timing studies, do not contradict inside the errors to its identification with G" star and 
' its proper motion. These errors are larger then ones in optical measurements, but are smaller then 

' corresponding errors in X-ray and 7-ray data. Estimations of the gamma-ray pulsar coordinates and 

its proper motion could be obtained independently on its optical or radio component, and are available 
in their absence. 

Q . 1 Introduction 

. The discovery of Geminga as a "true" pulsar, but without visible radioemission, gave ad- 
ditional evidence to the idea that hard gamma-ray emission {E > 30 MeV) is an inherent 
property of pulsar radiation. Before this only the Vela pulsar, the strongest hard gamma- 
^ ' ray source and the young Crab pulsar gave hints of this possibility. Observations on the 



■ sky, made by EGRET on CGRO ||28| |^ have shown that only young pulsars with age not 



exceeding several tens thousand years and, possibly, millisecond pulsars |^ give observable 



flux in gamma radiation. It is not yet clear, what is the mechanism of this gamma radiation 
and whether it has a threshold character, or if there is a gradual decrease of hard gamma 
ray flux with an age. 

The existence of the Geminga pulsar indicates, that there could be other gamma ray 
pulsars with no radioemission, which are exhibited in EGRET observations as ordinary 
point-like sources, see i.g. |^, |1T|. 

Determination of pulsations in a hard gamma ray source is a very difficult problem, 
connected with rareness of arriving quanta 6t ^ P, and small total number of quanta. 
When the value of the period is known from other observations (radio or X-ray), timing 
analysis gives the possibility to reproduce this periodicity also in gamma region p5 



When there is no information about the period, it could, in principle, be found from pure 
gamma data , but this could take enormous amount of computer time and has not been 
in full realized in practice. 
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A position on the sky of pulsars with no radioemission cannot be estabhshed precisely; the 
best position obtained from X-ray observations are between 3" (Einstein) and 5" (ROSAT) 
for 90% level ^ |2[. Optical identification of Geminga with very faint > 251*5 object have 
been done in 0] and later measurements of its proper motion [^] and parallax |]13[ can be 
considered as an evidence of reality of this identification. 

Timing of Geminga in hard gamma region based on COS-B 0, and EGRET |25 



data gave anomalously high braking index = z/i>/z>^ ~ 10 — 30, corresponding to very high 
second derivative i). While there is a possibility, that it is connected with poor precision of 
i) determination, it is worth to investigate other explanations. It was suggested in [0], that 
high value of n (i)) results from errors in its coordinates, leading to incorrect barycenter 
reduction procedure, which spoils the timing procedure. This problem is well known for 
pulsar timing, where the error in coordinates give a one year periodical deviations from 
the smooth curve in the pulse arriving time, what permits to improve pulsar position to 
amazing precision of the order and even better than VLBI observations [0, [27[] . 

Here we describe a method for investigation of timing of gamma pulsars, represented by 
periodical objects with rare pulses, which gives possibility to determine 7 parameters of a 
gamma pulsar: frequency z/, its two derivatives u and i), angular coordinates a and 6 of 
the source, absolute value v and direction of a velocity of a proper motion, characterized 
by an angle 6. When registered pulses are rare, so that their time separation 6t is much 
larger then the period P, the method of investigation is quite different from the one, used 
for radiopulsars. On the artificial sample of data, which properties simulate Geminga, but 
in contrary, have very narrow light curve (5-function), no systematic errors and no false 
quanta, we have managed to determine coordinates and proper motion parameters with 
very high precision. This precision is decreasing when we go to pulses with a finite width, 
in presence of background and systematic errors. 

Application of this method to real data sample of Geminga from COS-B and EGRET 
was not so successful, because of smooth light curve, "nonperfectness" of data, possible 
glitches. We present here results for most probable position and proper motion charac- 
teristics, determined by using of periodicity criteria, which are not in contradiction with 
more precise optical data, and has better coordinate precision then 7 - ray or X - ray data. 
Correlation properties of timing criteria, used for sources with rare pulses are investigated, 
before applying them to a timing procedure. 



2 Barycentric corrections: account for angular coordinates cor- 
rection and proper motion 

For timing analysis all data must be presented in the same coordinate system, which as a 
rule is connected with a barycenter of the Solar system. Consider first a situation, when 
angular coordinates of a source a and 6 are known exactly. On Fig. 4 xyz is a coordinates 
system, that remains at rest with respect to distant stars. The point O is the barycenter of 
the Solar system. The space probe is in the point S with Cartesian coordinates {xo,yo,zo). 
A direction to a source is defined by a straight line with coordinate angles a and 6, where a 
angle is counted in xy plane counterclockwise from the positive direction of a;-axis. An angle 
6 is countered in a plain, perpendicular to xy plane. SB is a perpendicular from the space 
probe (point 5*) onto a line from the barycenter to a source. BC is a perpendicular from the 
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point B to the xy plane, and CD is a line parallel to OB; points C and D belong to the xy 
plane. The unit vector pointed from the barycenter to a source has Cartesian coordinates 
(cos a cos 6, sin a cos 6, sin 6) . The scalar product of this vector and the radius- vector of 
the space probe is the length of the segment OB. 

OB = xq cos 6 cos a + yo cos 6 sin a + zq sin S 

A time interval during which light flies the length OB is a barycenter correction, if a source 
is far enough and SB is a part of a fiat wave front. This time interval AT, which must 
be added to the moment of each event in the point S to obtain a corresponding barycenter 
moment, is determined as 

AT = - (a;o cos 6 cos a + yo cos 6 sin a + zq sin 6) (1) 

In a common choice a and 6 coincide with right ascension and declination, when xy is an 
Earth's equatorial plain for some fixed epoch and x-axis points to spring equinox at the 
same epoch. Procedure of calculation of AT with precise account of Earth and satellite 



motion is described in 24 



Suppose that coordinates a and S are known not exactly with corresponding errors da 
and d6. Then for small errors we may find from (1) corresponding barycenter corrections in 
linear approximation 

dAT dAT 
5T = — — da H — 7:T—d6 = 
^oa do 

= -(^— Xq cos 6 sin a + yo cos 6 cos a^ da + (2) 
+ - Xq sin 6 cos a — yo sin 6 sin a + zq cos dS 

Consider for simplicity a case when a space probe orbit around the Sun is circular and, 
consequently, its angular velocity u; is a constant. Then xq = Rcosut, yo = Rsinut, ^0 = 
and 

R R 

AT = —(cos cut cos 5 cos a + sin cut cos 5 sin a) = — cos5cos(u;t — a) (3) 
c c 

where R is a. radius of its orbit. The correction (2) then is reduced to 



ST 



R 



cos 6 sin{u!t — a) da — sin 6 cos{ujt — a) d6 . (4) 



Assume that only first and second derivatives of u are essential, and a frequency of the 
signal may be represented by 

u = uo + uot + -j-f, (5) 

where index "0" is referred to the epoch t = 0. Define a current arrival time of photons 
from the source, measured on the satellite, as t. When the source coordinates are known 
exactly, barycenter arrival time t is found as 

t = t + AT, (6) 
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and having barycenter arrival times it is possible to find vq, z>o and i>o using the criteria from 
previous sections. 

When source coordinates are not known exactly and their possible errors are da and dS, 
the error in barycenter correction is determined by (4). In order to estimate an input of 
these errors on timing characteristics let us compare phases of the arriving signal calculated 
from measurements (f)' (with errors) and in true barycenter time 0, so that 

(f)' = f u'dt', (j)= f' udt (7) 
JO Jo 

Here t' is the time calculated from (6), and u' is a frequency found after barycenter correc- 
tions, containing errors. Times t and f correspond to the same event, so we may rewrite 0' 
in true barycenter coordinates as 

Using (4) and (5) in (8) we obtain after integration and account of (4), (5) 

0' = const + / iydt + v5T - [ 5T{vq + i>oi)di 
Jo Jo 

= const + vdt +[v-^5T- (z>o + i)ot) STdt. (9) 
Here and farther relations 

f r d^T f 

ST^-u^ST, J {J STdt)dt^-^, ^^-cu^JSTdt (10) 

are used, and it follows from (4) 

/R 
5Tdt — [cos(5cos(a;i — a)da + smSsm{u}t — a)d5]. (11) 
CUJ 

Differentiating (9) we obtain an input of the angular coordinate errors into the values of 
frequency and its derivatives 

dt \ dt J ' 
i/ 1 + — - - vuj dT, (12) 



dt^ \ dt 

dt3-^+(^---') dt 



u — -— — u + [u — uuj )— 21/00 oT, 



where 



z> = z>o + i>ot, i> = i^o, (13) 
and current values i/', v' and v' are connected with corresponding values at i = as 
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= ^0 + i'ot + i)'^- 

K = v'-v't^v'\ (14) 
z>^ = v' - v't 
v', = l)' 

The detailed variant of previous calculations can also be found in [^. 
In presence of a proper motion of the source the errors da and d6 change linearly in first 
approximation as 

da = dao + Vat, d6 = d6o + vgt. (15) 

It leads to farther complication of the formula (9)- (12). Note that in simulations we deal not 
with these formulae, but directly with arrival times of photons tj, t'^ and ti. In the problem 
of timing of radiopulsars the precision of observational data is very high, so appearance of 
the periodical 1 year component gives a direct indication to the errors in angular coordinates 



of the pulsar, possibility to improve them |2^, |T9|, and to determine a proper motion. In 



periodic sources with rare pulses a quality of data is much worse and other methods, based 
on above mentioned criteria must be used. 



3 Mathematical simulation 

For checking a possibility to use criteria considered above for determination of corrections to 
the angular coordinates and proper motion, in addition to frequency and its two derivatives, 
artificial sample of data was produced. A pulse shape was taken as (5-function with a 
frequency of the signal changing in time according to (5), what corresponds to a phase 
dependence 

(j) = (f)o + uot + i>o— + i)o — . (16) 

/ 

We need to find time moments, corresponding to phase values (pi = 27ri. Two sets of input 
parameters were considered. The time t = is related to a point of the orbit, where a = 0. 

(i) 00 = 0, z/Q = 4s-\ z>o = -2 ■ 10"®s-^ uq = 3 ■ 10'^^ s'^, 

= 6.060171 ■ 10"^-^ (17) 

(ii) 00 = 0, 1^0 = 4:S-\ z>o = -2 ■ 10^l3s^^ = 3 ■ 10-2S~^ 

CO = 1.991063802- 10"^ s"^ 

The source parameters are chosen to satisfy a relation vqVq/vq = 3, supposed to be valid for 
ejecting pulsars [^. Because of low reliability of i) detecting in gamma observations, some 
authors set it equal to zero. The modehng year duration (27r/u;) is equal to 12 days in 
the first case; and is a true value of 365.2422 days in the second, when the moment t = 
corresponds to 21 March. Note that in the second case the values of uq and i>o are chosen 
very close to that of Geminga [ED . 
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One possible way to find ti is to use Burmann-Lagrange expression, which links the Taylor 
coefficients of direct and inverse functions. We have used instead a procedure, valid for a 
general law of a phase dependence 0(t), based on a Taylor expansion formula 

To find the coefficients An = n> 1, use an evident equality 

• dt 



1. (19) 



That gives 



dt 1 

^1 = 71 =T- (20) 
Differentiating (19) (n — 1) times over t we obtain a relation, linear to An, what permits to 

[m] [n] 

express An as a function of Am, m < n — 1, and 0,0 . As an example, after 5 differentiation 

we get 

[6] [5] [4] [4] 

Ai0o + 6^20000 + 15^20000 + 10^20^ + 15^30000 + 60^3000000 

+ 15^301^ + 2OA40ii0o + 45^40^0^ + 15^50^00 + ^600 = 0, (21) 

[m] 

where index "0" indicates time t = 0. For (16) with 0o = at m > 4 we have 0o = I'o, 
00 = i^o, 00 = i^Q and get from (21) an equation for Aq 

10^2^0 + 60A3Z/oz>oi>o + 15A3Z>^ + 20v44Z/^^o + ^^A^u^u^ + l5A5U^uo + Aqu^ = (22) 
The first criterion Ki of periodicity [|l^, |16|, |l^, [l^ was used to investigate periodicity 



properties of series of pulses. For a purpose of testing short intervals of "observation" were 
taken in a different parts of the year. The error in coordinates was taken equal to 5 arc 
seconds in absolute value ^ dcP' + db'^, but the deviations from the initial point were taken in 
eight different directions, separated by 45°. The values of v'^, z>q, v'^ that have been detected 
by Kx criterion coincide in both cases with a very high precision with theoretical ones from 
(12)-(14), see Table 1. It may be seen from Table 1 a strong inffuence of the errors on the 
determination of v' by using a criterion. While the error in i)' is almost linearly proportional 
to the error in V dd^ + db'^, it is evident that at an angular error larger then 10"'^ arc sec 
direct determination of v' by criterion becomes impossible. This may be a reason for a high 
breaking index of Geminga pO|, 0. 

Let us now formulate a problem of timing of a gamma pulsar, which gives a possibility 
for a search of its timing properties together with angular coordinates and a proper motion. 
Assume that a gamma pulsar simulated by computer is emitting signals with a frequency 
changing according to (5), satisfying condition vt^vt^jv^ = 3. 

Let the signal registered on the probe is reduced to the barycenter time, using the source 
coordinates ao and 6q (base point), which contain errors da^ and d6o respectively. We 
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suspect also a proper motion of the source defined by following parameters: at the moment 
to = the source has coordinates ao + dao and 60 + d6o and a velocity modulus 9 = const. 
A velocity direction is defined by an angle 6, which is counted clockwise from the positive 
a-axis direction. The current source coordinates are consequently 

a = ao + (iao + ^^<^ cos6', 6 = 60 + d6o + v^tsinO, (23) 

Thus, we have seven parameters that are needed to be found selfconsistently 
uq - frequency of the source signal 
z>o - first derivative of the frequency 
i>o - second derivative of the frequency 

dao ^ shift in right ascension from the base point (ao, Sq) at to epoch 
d6o - shift in declination from the base point (ao, 60) at to epoch 

- proper angular velocity of the source in celestial coordinates 
6 - direction of velocity , counted clockwise from the positive a-axis 
direction 



with additional restriction: z/o^'o/^'o ~ 3 following from the model of the pulsar radiation p2 

In data simulation we fix parameters z/q, z>o, i^o as (ii) in (17), find true a and 6 from 
(23) with dao = 2", d6o = 3", = 0'.'2 per year, 9 = 60°, and create the simulated data, 
i.e. the sequence of time moments of pulses. Time moments (true barycenter) found for a 
source from (19)-(22) are then recalculated for a probe using correct coordinates and (23). 
Now to a set of time moments "registered" by a probe from the source with subscribed 
coordinates ao and 60, containing errors, we apply the algorithm for searching the periodic 
signal to extract all the seven parameters from the simulated data set. Namely, we consider 
a number of different sets of 7 mentioned parameters. For each set we evaluate supposed 
errors, introduced in the data due to errors in a position and in a proper motion of the 
object. After that we subtract these supposed errors from the data. Then, assuming the 
data free of errors the periodicity criterion value was calculated. Remind, that if we assume 
the data free of errors, it means that the phases of the pulses must obey the simple relation 
(16). The criterion reaches its absolute maximum only for an exact set of parameters. There 
appear a number of local (or false) maxima. Fig. 5 demonstrates a typical structure of the 
criterion depending on two parameters: u and u, when other five ones are fixed (see also 
[^). Because the height of a "false" maxima is close to "true" one, it is very difficult to 



separate the absolute maximum among a series of local ones. 

We are looking for an absolute maximum, using a grid in 4-dimensional space 
{dao,d6o,v^,6), that was defined by the following way: dao varies from —5" step O'.'S to 
5" , d6o varies from —5" step O'.'S to 5" , varies from 0"/year step O'.'05/year to O'.'5/year 
and 6 varies from 0° step 10° to 350° . Three other parameters (z/o, z>o, i)o ) were detected 
jointly for each point of above grid, using the grid 6x6x6 with 12 times consecutively 
diminishing steps for all three axes. The maximum, detected at the previous step was placed 
to the center of the grid and all the scale multiplied by the factor of 0.4. This procedure was 
repeated for 12 times, so that the last grid steps are 0.4^^ ^ 1.7 ■ 10~^ times as small as the 
first ones. All seven parameters chosen for modeling were found with a precision, limited 
only by computational grid connected with a power of the computer, using criterion Ki. So, 
for a clean set of data the proposed procedure of searching is working effectively. Situation 
is becoming much more controversial when we apply it to real data existing to the moment. 
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4 Application to Geminga 



4.1 Analysis of COS-B data 

COS-B mission had been operated since August, 1975 till April 1982, and had observed 
Geminga in five shifts p^, ^. Their numbers are 00, 14, 39, 54 and 64. Some measurements 
of Geminga position, reduced to the epoch 1950.0, are summarized in Table 2 and are plotted 
in Fig. 6. 

It was obtained in |g from COS-B data v = 4.217 Hz, i> = -1.952 ■ lO-^^n^.g-^ and a 
large value of i) = (28 ± 16) ■ 10^^^ Hz-s~^, corresponding to a braking index n = 31 ± 18. 
Barycenter corrections have been done for standing Geminga with coordinates No. 3,8 in 
Table 2. Quanta selection used in our analysis has been done by two different ways: 

1) all the quanta in the circle of r = 5° around Geminga position; the interval 54 

was excluded because of low reliability; it is total of 1505 quanta. 

2) all the quanta of the energy E > 50 MeV laying in the circle r = 12.5 ■ E~^'^^, 

where E is measured in MeV and r in degrees |T^; total of 1883 quanta. 



The second selection is close to that used in [^]. A problem of a quanta selection criterion is 
a very delicate one because it is practically impossible for a single quantum to decide, was 
it really radiated by Geminga or belongs to a background. Both mentioned selections are 
noisy, but the first one is worse. 

Geminga was considered as a moving object. The base coordinates, that are used for 
initial barycenter reduction are the position measured by Einstein satellite in 1981 (see 
Table 2). Geminga motion was defined by its velocity, direction and initial position at the 



epoch 1979, March, 14.0 ||20[. The obtained barycenter time moments for each quantum 
was additionally reduced to the barycenter, using expression (2). A motion of the probe 
defined by xo(t), yo(t), zo{t) was taken from databases of COS B or EGRET. The object 
coordinates were calculated by this procedure separately for each quantum according to a 
supposed object motion. 

The results, with using the second selection from the mentioned above, are not very 
certain. The criterion appear to exhibit a gently sloping maximum at the following model 
parameters: velocity = 0'.'2 — 0'.'3 per year, direction 6 = 40° — 60° and the initial 
coordinate offsets dao and d6o at the mentioned epoch are —2" for both a- and 5-axes, but 
the uncertainty here is high and may reach 2" for both coordinates. As to the periodicity 



parameters, they are in a good agreement with except the second derivative uq, which 
is a bit smaller but lays within the error box of standing Geminga. The motion of Geminga, 
obtained in our investigation does not contradict to the motion of G" star . 

Unfortunately, this solution is not a unique one, and there are a number of other maxima 
of approximately the same height. When we search for a global maximum in 7-dimensional 
space it is extra difficult to detect "the main maximum" among a series of other local 
maxima. For example, there is an accessory maximum at = O'.'l — 0'.'2 per year, 6 = 
340° — 360° and very badly detected initial offsets (it can only be said that they both are 
negative). The periodicity parameters here are approximately the same as above, but uq 
appears to be negative. 
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4.2 Analysis of EGRET data 

EGRET experiment is operating since April 1988. There are 9 periods of observations 
where Geminga was not far from the center of a view field (less than 30°, the standard 
requirement). The following sessions was used for data investigations: 2, 3, 4, 5, 10, 21, 
2130, 2210, 3100. The standard procedure for the barycenter correction was used P^, P3 



but we have used coordinates of the object from [Q, line 8 in Table 2, different from those, 
indicated in EGRET data base. We have used the same Geminga coordinates for both 
satellites, COS-B and EGRET. They are the best fit Einstein position, but in the last case 
we were to reduce them to the epoch 2000.0, because it is used in an appropriate barycenter 
reduction routines. Reduced coordinates are the following: 

^2000 = 98°28'30'.'90 = 98?47525 ^2000 = +17°46'11'.'6 = +17?76989 (24) 

And in EGRET data base the coordinates are: 

^2000 = 98?48 = 98°28'48" ^2000 = +17?77 = +17°46'12" 



Other authors |jT^, ^ use the coordinates close to (24), they differ less than 1" from the 
center of Einstein error box (lines 3, 8 in the Table 2). 

We have used a number of techniques for quanta selection and have compared the results. 
The selections used are the following: 

1) all the quanta of the energy E > 70 MeV laying in the circle r = 5.85 ■ 

{E /100)~^'^^^, where E is measured in MeV and r in degrees it is so- 
called the standard selection; total of 6751 quanta. 

2) all the quanta of the energy E > 1500 MeV in the circle r = 2°; total of 365 

quanta. 

3) all the quanta of the energy E > 2000 MeV in the circle r = 2°; total of 223 

quanta. 

4) all the quanta of the energy E > 2000 MeV in the circle r = 0?5 round the 

Geminga position; total of 100 quanta. 

5) all the quanta of the energy E > 3000 MeV in the circle r = 0?5 round the 

Geminga position; total of only 51 quanta. 

The criterion value for unmoving Geminga, resting in Einstein's HRI position is Ki = 
0.0572 and the parameters of periodicity are u'^ = 4.21775012925 Hz, z>^ = -1.95312 ■ lO'^^ 
Hz-s-\ u'q = (20 ± 12) ■ 10-26 Hz-s-2. An appropriate light curve is shown in Fig. 7. The 
results of search in 7-dimensional space for different selections are: 

1) There is a gently sloping maximum in criterion value at the following parameters: z/q = 
4.2177501295 Hz, z>o = -1.9532 • lO^^^ Hz-s-\ i)o = (22 ± 13) ■ lO^^e Hz-s^^ f< = 
0'.'3 — 0'.'4 per year, 6 = 40° — 60°, dao and d6o are defined with very low precision 
and both lay in the interval —2" to 0". The criterion value Ki = 0.05732. The second 
derivative i) here is of rather large value, so that the braking index is approximately 
equal to 25. 
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2) There is a gentle maximum in criterion value at the following parameters: Uq = 

4.2177501273 Hz, z>o = -1.952711 ■ IQ-^^ Rz-s'^ />o ~ (0 ± 10) ■ IQ-^^ Hz-s'^ = 0'.'3 
per year or more, 6 = 40° — 60°, dao = 0" — (—2"), d6o = 0" — (—1") and criterion 
value Ki = 0.22357. The second derivative here is very close to zero and the braking 
index is, respectively, also low and could be close to its theoretical value. 

3) There is a gentle maximum in criterion value at approximately the following parameters: 

z/o = 4.2177501261 Hz, z>o = -1.952611 ■ lO'^^ Hz-s-\ i)o = (-2.4 ± 10) ■ 10-^6 Hz-s-2, 
= 0'.'3 - 0'.'4 per year, ^ = 40° - 80°, dao = 0" - (-2"),rf5o = -1" - (-3") and 
the criterion value Ki = 0.25535. It was impossible to determine the parameters with 
higher precision. 

4) There is a gentle maximum in criterion value at the following parameters: uq = 

4.2177501344 Hz, z>o = -1.954477 ■ lO-^^ Hz-s-\ i)o = (68 ± 40) • lO'^e Hz-s'^, 
= 0".3 — 0"A per year or more, 6 = 220° — 260°, dao and dSo are negative and 
the criterion value Ki = 0.3072. Because of a very poor statistics this selection as well 
as the following one can be considered as a test only. The line of motion here is the 
same as in the previous cases, but the direction is opposite . 

5) There is a relatively good maximum despite of a very poor statistics at the following 

parameters: z/q = 4.2177501220 Hz, z>o = -1.95282-10-13 Hz-s~\ uq = (12 ±10) • lO'^^ 

Hz-s-2, = 0'.'2-0'.'3 per year, 9 = 20°-60°, dao = 0" 2", d6o ^ 0" and the criterion 

value Ki = 0.44441. The reliability of this result is not high, but large criterion value 
indicates that we may deal with real motion of the object. 

4.3 Analysis of combined COS-B and EGRET data 

We have combined COS-B and EGRET data with the following selection criteria: 
COS-B: ^ > 50 MeV and the standard conditions for r: r = 12.5 ■ E'^-^^ [|l|] 
EGRET: E > 70 MeV and the standard conditions for r: r = 5.85 • {E/100)~°-^^^ [g9| 



There are total of 1883 + 6751 = 8634 quanta. 

The combined series is not self-contradicting. There is a gentle maximum in crite- 
rion value at the following parameters: = 4.21775012323 ± 0.00000000025 Hz, z>o = 
(-1.952554 ± 0.000025) ■ lO'^^ s'^ i)o = (-2.5 ± 10) ■ lO'^^ Hz s'^, = 0".5 - 0'.'6 per 
year, ^ = 55° - 65°, dao = -2" - (-4"), dSo = 1" - 2" and the criterion value Ki = 0.04937. 

Note that la errors above were obtained by approximate estimations. 



5 Discussion 



According to our investigations the coordinates and motion of Geminga obtained from 
timing of gamma pulsar is in a satisfactory agreement with the motion of G" star when 
separately COS-B or EGRET data are used. Parameters following from the combined data 
set are in much worse agreement. There are two possible reasons of it. First there could 
exist a systematic error between the data of two probes; and second, period could behave 
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nonmonotonously between 1982 and 1988 years, and the period jump (pulsar glitch) of the 
order of ^ > 10~^° could already spoil the parameters obtained by criteria. 

The value of the second derivative i)o does not coincide with the theoretical one, however 
it is lower than in the previous investigations, but there are weighty reasons to explain 
this phenomenon. The second derivative is a very sensitive variable, and even O'.'OOl error 
in angular coordinates changes uq by the value, comparable with the result (see Table 1). 
If there was a jump in pulsar period, it may also cause the incorrect value of a variable. 
Possibility to improve angular resolution by timing is strongly limited by small number of 
quanta and existence of considerable background. 

As a result of application to Geminga of the developed method of timing of gamma 
pulsars we have obtained that determination of true value of i)o is possible only at very high 
precision (better then O'.'OOl) of angular localization. At good statistics of gamma pulsars 
corresponding improvements would become possible from timing analysis. New types of 
gamma ray telescopes based on very wide aperture (> 2.5 vr steradian) and higher threshold 



of a few hundred MeV 21] would permit to get higher angular resolution (~ 1 arc min) 



reducing influence of a background, get ~ 100 better statistics due to continuous monitoring 
of larger part of the sky in this region. 

For the existing data of Geminga from COS B and EGRET it was obtained, using only 
gamma-ray data, that criterion value reaches its maximum at nonzero value of a proper 
motion. The coordinates of the source were confirmed with precision ~ 2"what is better 
then follows from X ray data, but, of course, is worse then the precision obtained in the 
optical observations. 

Observations of radio pulsars have shown, that their optical and X-ray luminosity is de- 
creasing with time much more rapid, that radio and hard gamma radiation. So at increasing 
sensitivity we expect a discovery of tens of new gamma-ray pulsars similar to Geminga, may 
be without X-ray and optical counterparts. For such objects method of timing of gamma- 
ray pulsars developed above would be a main and may be a single means of investigation of 
such sources by data processing. 
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Table 1. 

Properties of a source with angular coordinates and 

timing characteristics close to Geminga from (17), case (ii), 

calculated using criterion (1). Barycenter corrections 

have been done using angular coordinates with indicated errors. 



Source coordinates: a = 6''30™00^ 5 = +17°30'00" 


Montn, 

day 


Errors 


Theoretical values 


da, 


do. 


—1 

I/O, s 


•~ 1 A— 13 —2 

Uq, 10 s 


■•' 1 n— 26 —3 


May, 20 


5.00 


0.00 


3.99999999748 


-1.9738 


-28600 


May, 20 


2.24 


2.24 


3.99999999744 


-1.9879 


-9620 


May, 20 


0.00 


5.00 


3.99999999681 


-1.9992 


7140 


May, 20 


-2.24 


2.24 


3.99999999970 


-2.0114 


16000 


May, 20 


-5.00 


0.00 


4.00000000252 


-2.0262 


28600 


May, 20 


-2.24 


-2.24 


4.00000000256 


-2.0121 


9620 


May, 20 


0.00 


-5.00 


4.00000000319 


-2.0008 


-7140 


May, 20 


2.24 


-2.24 


4.00000000030 


-1.9989 


-16000 


July, 10 


5.00 


0.00 


3.99999999621 


-1.9695 


-35700 


Aug., 20 


5.00 


0.00 


4.00000000544 


-1.9852 


-22300 


Aug., 20 


0.00 


-5.00 


4.00000001020 


-2.0155 


9110 


Nov., 20 


5.00 


0.00 


4.00000008067 


-2.0725 


29150 


Nov., 20 


-2.24 


2.24 


3.99999996211 


-1.9631 


-16200 
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Table 2. Geminga position from different measurements 



Geminga position, epoch 1950 


No. 


a 


b 


Error 


Comments 


1 


97°44'43'.7 


+17°48'27'.'5 


12" 


ROSAT PSPC, 1991, Sep. 19-21, 


2 


97°44'51'.7 


+17°48'36'.'0 


^ 5" 


ROSAT HRI, 1991, Mar. 19, 


3 


97°44'47'.'2 


+ 17°48'33'.'0 


3"2 


Einstein, 1981, Mar. 18, [§ 


4 


97°44'45'.'9 


+17°48'32'.7 


0'.'46 


G" star, 1984, § 


5 


97°44'45'.'9 


+ 17°48'32'.'6 


0"5 


G" star, 1986, Feb. 3, § 


6 


97°44'46'.'5 


+ 17°48'33'.'0 


0"68 


G" star, 1987, § 


7 


97°44'47'.'2 


+ 17°48'33'.'6 


0"16 


G" star, 1992, § 


8 


97°44'47'.'2 


+17°48'33'.'0 


3"0 


Einstein, 1981, Mar. 18 [§ 
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Pulses 
400 J 




Epoch 1979, March, 14.0 
i>o = 4.21775012925465 Hz 
z>o = -1.95316166 • 10"^^ Hz s" 
i>o = 19.66 • 10-26 Hz s-2 
n — vqi/q/vq — 21.73 



Figure 4: Geminga 40-biii light curve from ERGET data, for standing Geminga with cooordinates No. 8 
from Table 2 and standard quanta selection 
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Figure captions 



Figure 1: On the barycentric correction. 
Figure 2: On the absolute maximum structure. 
Figure 3: On the Geminga position. 

Figure 4: Geminga 40-bin hght curve from ERGET data, for standing Geminga with 
cooordinates No. 8 from Table 2 and standard quanta selection. 
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Figure 1: On the barycentric correction 
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Figure 2: On the absolute maximum structure 
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Epoch 1950 



Figure 3: On the Geminga position 
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